Spin Waves in Striped Phases 
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In many antiferromagnetic, quasi-two-dimensional materials, doping with holes leads to "stripe" 
phases, in which the holes congregate along antiphase domain walls in the otherwise antiferromag- 
netic texture. Using a suitably parametrized two-dimensional Heisenberg model on a square lattice, 
we study the spin wave spectra of well-ordered spin stripes, comparing bond-centered antiphase 
domain walls to site-centered antiphase domain walls for a range of spacings between the stripes 
and for stripes both aligned with the lattice ("vertical") and oriented along the diagonals of the 
lattice ("diagonal"). Our results establish that there are qualitative differences between the ex- 
pected neutron scattering responses for the bond-centered and site-centered cases. In particular, 
bond-centered stripes of odd spacing generically exhibit more elastic peaks than their site-centered 
counterparts. For inelastic scattering, we find that bond-centered stripes produce more spin wave 
bands than site-centered stripes of the same spacing and that bond-centered stripes produce rather 
isotropic low energy spin wave cones for a large range of parameters, despite local microscopic 
anisotropy. We find that extra scattering intensity due to the crossing of spin wave modes (which 
may be linked to the "resonance peak" in the cuprates) is more likely for diagonal stripes, whether 
site- or bond-centered, whereas spin wave bands generically repel, rather than cross, when stripes 
are vertical. 



I. INTRODUCTION 



Many doped strongly correlated materials exhibit ev- 
idence for an emergent length scale in the form of 
"stripes", i.e., regular antpihase domain walls in an 
otherwise antiferromagnetic texture. The strongest ev- 
idence for striped structures in nickelate perovskites 
and some related cuprates has come from neutron 



scattering, 



1.2.3.4.5 



which is capable of detecting the spin 



texture directly through diffraction. Since several theo- 
ries of high temperature superconductivity make contact 
with such structureSfSi^iSiiSiiiii^ it is important to im- 
prove our microscopic picture of them. In particular, it 
is not yet known from experiment whether the antiphase 
domain walls sit primarily on nickel (copper) sites, or 
rather sit primarily on oxygen sites. 

When undoped, the nickel-oxygen (and copper- 
oxygen) planes in these materials are antiferromagnetic, 
with spin moments localized on the Ni (Cu) sites, as evi- 
denced by a peak in elastic neutron scattering at (ir, ir)i& 
Upon hole doping, this peak is observed to split into four 
(or perhaps twoi^ii 5 .) "incommensurate peaks" indi- 
cating an extra modulation on top of the antiferromag- 
netic wavelength. For the case of colinear spins, this is 
consistent with the formation of periodic antiphase do- 
main walls in the antiferromagnetic texture (i.e. stripes). 

On a two-dimensional square lattice, these domains 
consist of a strip of antiferromagnet with spin up on, 
say, the "A" sublattice, separated by a domain wall from 
a strip of antiferromagnet with spin up on the "B" sub- 
lattice, and so on, as in Fig. ^ The figures necessar- 
ily depict a certain width for each antiphase domain 
wall, but the widths are not known and are in real- 
ity likely less sharp than shown in the figure. In both 
cases, neighboring antiferromagnetic patches have spin 



up on opposite sublattices, which washes out any sig- 
nal at the antiferromagnetic peak (jr, it). Rather, satel- 
lite peaks are observed around (71", tt), at a distance de- 
termined by the spacing between domain walls. When 
the domain walls are site-centered, all couplings are an- 
tiferromagnetic, including couplings across the domain 
walls. Bond-centered domain walls, however, have some 
ferromagnetic couplings. 17 That is, bond-centered config- 
urations consist of antiferromagnetic patches which are 
ferromagnetically coupled across the domain wall. As 
shown in Fig. ^ we consider stripes aligned with the lat- 
tice direction (called "vertical stripes" ) or aligned along 
the lattice diagonals (called "diagonal stripes"). 



In this article we focus on the spin wave spectra and 
expected magnetic scattering intensities of bond-centered 
and site-centered stripe phases of various spacings and 
orientations. Other stripe phases are certainly possible, 
such as phases which mix site- and bond-centered do- 
main walls, or phases in which the spacing of the an- 
tiphase domain walls is not commensurate with the un- 
derlying lattice, or "dynamic" stripes^, which fluctuate 
in time, We will not consider these cases here, but fo- 
cus on well-ordered spin stripes which have purely site- 
or bond-centered domain walls. As we will show below, 
there are qualitative differences between the spin wave 
spectra of bond- and site-centered domain walls, indicat- 
ing that in some cases inelastic neutron scattering may 
be able to distinguish between the two. In addition, there 
is a difference in the number of peaks in the elastic spin 
structure factor for odd stripe spacings, indicating that 
elastic neutron scattering alone may be able to distin- 
guish as well. 
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FIG. 1: (Color online) (a) Site-centered vertical stripe pat- 
tern with p — 4 lattice constants between domain walls. In 
this configuration, exchange couplings J a > and Jb > 
are all antiferromagnetic. (b) Bond-centered vertical stripe 
pattern with spacing p — 4. The exchange coupling J a > 
is antiferromagnetic, while Jb < is ferromagnetic, (c) Di- 
agonal site-centered domain walls have coupling Jb > for 
next nearest neighbor spins coupled across the domain wall 
along the vectors (2,0) and (0,2), and coupling J c > diag- 
onally to "Manhattan" second neighbors across the domain 
walls along the vector (1,1), in units where the square lat- 
tice spacing a = 1. (d) Diagonal bond-centered domain walls 
have nearest neighbor ferromagnetic coupling Jb < across 
the domain wall. The size of each figure has been chosen for 
visual clarity. 



II. MODEL 



A. Vertical Stripes 

We consider first the case where stripes run parallel to 
the Ni-0 (Cu-O) bond direction; we call these "vertical" 
stripes. As illustrated in Fig.^ when the domain wall is 
centered on a lattice site, we may describe the system as 
having no net spin on the domain wall^ In this case, 
spins from the edges of neighboring antiferromagnetic 
patches are coupled across the domain wall antiferromag- 
netically, J r . r ' = Jb > with S r = on the domain wall, 
as illustrated in Fig. ^a). Within the antiferromagnetic 
patches, nearest neighbor spins are of course also antifer- 
romagnetically coupled, J r r ' = J a > 0. When, however, 
the domain wall is bond-centered — that is, situated be- 
tween two sites as in Fig. db) — spins from the edges of 
neighboring antiferromagnetic patches are ferromagneti- 
cally coupled, and we have J r r i = Jb < across the do- 
main wall. Nearest neighbor exchange couplings within 
each antiferromagnetic patch remain antiferromagnetic, 
Jr,r> = Ja > 0. We shall see that this ferromagnetic 
coupling Jb of spins across the domain wall leads to dis- 
tinctive features for the spin waves in the bond-centered 
case. 

We define the magnetic Bravais lattice as follows^ 
Let p denote the distance between domain walls. We will 
henceforth work in units where the square lattice spacing 
= 1. For p — odd, we choose the basis vectors Ai = 
(p, 0) and A 2 = (0, 2), and for p = even, we use Ai = 
(p, 1) and A2 = (0,2). For site-centered configurations, 
there are N = 2p sites within each unit cell which include 
2(p — 1) spins and 2 sites with no static spin component. 
For bond-centered domain walls, there are N = 2p spins 
in each unit cell. (See Fig. 0) 

We use the notation VSp and VBp to refer to vertical 
stripes of spacing p in a site(S')- or bond(£?)-centered 
configuration, respectively. For example, ^5*3 refers to 
a vertical site-centered configuration with spacing p = 3 
between domain walls. 



We consider static, ordered arrays of antiphase do- 
main walls in an otherwise antiferromagnetic texture. Al- 
though the domain walls collect chargei2i22i2ii2£ ; we ne- 
glect this charge component, as we are interested solely 
in the response of the spin degrees of freedom. We use a 
Heisenberg model on a two-dimensional square lattice: 

H= - >^r,r'S r S r /, (1) 

<r,r'> 

where < r, r' > runs over all spin sites, and the ex- 
change coupling is J r , r '- Within an antiferromagnetic 
patch, nearest neighbor couplings are antiferromagnetic 
with J r r i = J a > 0. Couplings across a domain wall 
depend upon the configuration and are enumerated be- 
low. All other couplings are neglected. When compar- 
ing to the nickel oxides (copper oxides), our lattice cor- 
responds to the nickel (copper) sites within the nickel- 
oxygen (copper-oxygen) planes. 



B. Diagonal Stripes 

For diagonal stripes, the antiphase domain walls are 
oriented along the (1,±1) direction in a square lattice 
(recall we have set the lattice spacing a — 1). For 
the same microscopic interaction strengths (deriving J r r i 
from, e.g., a Hubbard model), spins are more strongly 
coupled across the domain wall than in the vertical case. 
For example, with diagonal bond-centered stripes, each 
spin neighboring the domain wall interacts with two near- 
est neighbor (ferromagnetically coupled) spins across the 
domain wall, as shown in Fig. ^d). Contrast this with 
the vertical stripes of Fig.QJa) and (b), where each spin 
neighboring a domain wall interacts with only one spin 
across the domain wall. Diagonal site-centered stripes 
are even more strongly coupled, with two different types 
of interactions across the domain wall, one of which we 
label Jb because it connects spins along a bond direction 
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(connecting spins along the vectors (2, 0) and (0, 2) across 
the domain wall), and the other we label J c (connecting 
spins along the vector (1,1) across the domain wall), as 
shown in Fig.^c). 

For diagonal stripes, the magnetic Bravais lattice dif- 
fers from the vertical case. For p = odd spacing between 
domain walls, we choose the basis vectors Ai = (p, 0) and 
A2 = (— 1, 1), and for p — even, we use Ai = (2p, 0) and 
A 2 = (—1, 1). For site-centered configurations, when p is 
even there are N — 2p sites within each unit cell which 
includes 2(p — 1) spins and 2 sites with no static spin 
component, and when p is odd, there are N — p sites 
within each unit cell, which includes p—1 spins, and one 
empty site. For bond-centered domain walls, there are 
N = 2p spins in each unit cell when p is even, and there 
are N = p spins in the unit cell when p is odd. (See 
Fig. El) 

We use the notation DSp and DBp to refer to diagonal 
stripes of spacing p in a site (S)- or bond (-B)-centered 
configuration, respectively. 

III. SPIN WAVE THEORY 

The elementary excitations of ordered spin textures 
may be studied using the well-known technique of 
Holstein-Primakoff bosons. The same dispersion is ob- 
tained by quantizing the classical spin waves, and the 
methods are equivalent as S — > 00. We use each descrip- 
tion when convenient. As it is physically more transpar- 
ent, we review here the latter method^, discussing the 
former in Appendix IaI 

In the classical spin wave approach, each spin is treated 
as precessing in the effective field produced by its coupled 
neighbors, via the torque equations of a spin in a mag- 
netic fields The rate of change of the spin at position r 
is described by 

= Mr x HV/ , ( 2 ) 

where fi T and H^-' are respectively the corresponding 
magnetic moment and effective magnetic field at position 
r, defined by 

Mr = -.9MsS r 

n eff = ±Vj Sr , , (3) 
9VB ^ 

Within our model, Eqn. Q), the torque equations be- 
come 

J r _ _±(cy\^ j ,q z -Q z X^t ,<? y, i 

— r 2—l T < r r r 2-~l T ' T °r') 

r' r' 

^r ( Q z \ 7 ox qx \ y j qz \ 

—j— — —{ O r 2_^ J r,r'J r ' °r / , ^r'J r 'j 

r' r' 

f - »• 



where we have assumed large S and small oscillations, so 
that changes in S z can be neglected. We seek solutions 
of the form 

S* = Sf exp [t(k-r-wt)] 

S? = Sf exp [f (k ■ r - ut)] (5) 

where i labels spins within the unit cell, i.e. i = 
1, 2, • • • , N; N is the total number of spins in the unit 
cell; k = {k x ,k y ), and r = (r Xl r y ). Setting the determi- 
nant of the coefficients of Sf and Sf to zero yields the 
dispersion relations for the spin wave. 

We calculate the zero-temperature dynamic structure 
factor using Holstein-Primakoff bosons. 

S(k,co)=Yl E l</l^(k)|0>| 2 <5(c-c / ) (6) 

/ i=x,y,z 

Here |0 > is the magnon vacuum state and \f > denotes 
the final state of the spin system with excitation energy 
ujf. Since S z does not change the number of magnons, 
it leads to the elastic part of the structure factor. Single 
magnon excitations contribute to the inelastic response 
through S^k) and S w (k). 



IV. RESULTS FOR VERTICAL STRIPES 

We begin with our results for ordered, vertical stripe 
phases. We discuss magnon excitation energies as func- 
tions of momentum, the dynamic spin structure factors, 
the elastic response, the velocities of the acoustic bands, 
and analytic results for dispersion relations for small unit 
cell sizes. Fig. [21 shows schematic representations of ver- 
tical stripes that are site- and bond-centered, with both 
even and odd spacing. In this figure (in contrast to 
Fig. (JIJ) we have used the length of the arrow to rep- 
resent the net spin on a site. The net spin is expected 
to be smaller near domain walls (as it is always zero on 
a domain wall). Our zero frequency results incorporate 
this general form factor. For the finite u> results, we use 
a form factor with the same net spin on each occupied 
site. 



A. Elastic peak at (0, 7r) 

Elastic neutron scattering can in principle detect one 
important qualitative difference between bond- and site- 
centered stripes. For odd stripe spacings, both bond- 
and site-centered stripes have magnetic reciprocal lattice 
vectors at (0,7r). However, site-centered stripes are for- 
bidden from producing weight at (0,7r), whereas bond- 
centered stripes generically show weight at this point. 
This is related to the discrete Fourier transform of the 
spin structure. Taking advantage of the antiferromag- 
netic long range order in one direction and the finite 
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FIG. 2: (Color online) Vertical site- and bond-centered con- 
figurations, showing even and odd spacing. "S"' refers to 
site-centered configurations, and "B" refers to bond-centered 
configurations. The number label is the spacing p between 
domain walls. Dotted vertical lines mark antiphase domain 
walls. The solid boxes denote unit cells. The height of the 
arrows represents the net spin on a site, which is expected to 
peak between domain walls. 
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FIG. 3: (Color online) Schematic representation of vertical 
stripes with p =odd widths, indicating the pattern of the 
functions g(m) and f(n). Note that for the bond-centered 
case with odd stripe spacings, the function f(n) can have a 
net magnetization, producing elastic weight at the peak (0, 7r). 



spacing between stripes in the other, we can describe the 
spin structure in real space by a function 



S z (n,m) 



cos(irm) Aj 

3=0 

f(n)g(m) , 



(7) 
(8) 



where m is the discrete y coordinate parallel to the 
stripes, n is the discrete x coordinate perpendicular to 
them, and where f — p— 1 for p odd, with j' — 2p— 1 for p 
even. The functions f(n) and g(m) are shown schemat- 
ically in Fig. |21 The elastic scattering cross section is 
proportional to the Fourier transform of S z (n, m)»2& 

(%) K E el(fcmm+fc " Tl) (^( m ' n ))( 5Z (°' )) ( g ) 

\ / el ?ri n 



m j=0 u 



n i( — - ) r 

e v v ' 



= N m {S km ^ + 5 km _ w ) ]T A 3 H Nn6 k n ,-^f- ■ 
3=0 n 

We emphasize that this expression allows for any form 
factor and is not restricted to configurations where each 



occupied site has a full quantum of spin. In the case 
where each occupied site has the same net spin, the ra- 
tio of intensity at the main peaks (tt ± 7r/p, w) to that 
at (0, 7r) is 2 in the VB3 case, and 2.6 in the VB5 case. 
Site-centered stripes always have Aq — 0, while ylo is 
generically nonzero for bond-centered stripes (although 
it can be fine-tuned to zero). A finite j = term pro- 
duces elastic weight at (0,tt). This can be understood 
heuristically from considering the function f(n), shown 
schematically for the VB5 case in Fig. [3J Odd-spacing 
bond-centered stripes generically have a net magnetiza- 
tion in the function f(n), while symmetry forbids this for 
site-centered stripes. 



B. Analytic Results for small p 

For small p, which corresponds to small unit cell sizes, 
we can obtain analtic results for the dispersion relation 
of the acoustic mode. For the case VS3, we find 



/ LOVS3 



4(A + 1) + C-4(A + 1)L> (10) 
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where 



X= 0. 4 



1= 1 



I = 2. 5 



A = A 



(ii) 



C 



2\f(3k x ) + 8f(k y )-4f 2 (k y ) 



D = \l-f(k y )\jl- 



2\f(3k x ) 



» (A + l) 2 
and the function / is defined as 

/(x) = 1 — cos(a;) . 



(12) 



The acoustic spin wave velocity parallel to the stripe 
direction (k| \y) may be obtained by setting k x = above, 
and taking k y « 1. In this case, f(k x ) — 0, f(k y ) — > 



i ;„2 



fc y , and wys2 — ► where 



1 



(13) 



and vaf — 2^/2 J a S is the velocity of the pure antiferro- 
magnet with coupling J a and no antiphase domain walls. 
The spin wave velocity perpendicular to the stripe direc- 
tion may be similarly obtained: 



Vj_ 



3V2 /A(A + 3) 



A + l 



VAF- 



(14) 



For A >> 1, these approach vj_ — > (3/\/2)\/A vaf and 

V|| -> (1/2)VAuaf- 

For the case VB2, the problem reduces to diagonalizing 
a 4 x 4 matrix, with the result 

(t?) = 2(A " + 3A + 2) + ^ ~ 2 ^ (A2 + 3A + 2 ) 2 + s 

(15) 

where 

A = 2/(2fc„) (16) 
1 



B 



X 2 f(Ak x ) - Af(ky) - 4(A^ + 3X)f(2k x ) 



-4/(fc„)(l - f(k y ) + (A 2 + 3A)(1 - f(2k x ))) 
The spin wave velocities in the case VB2 are 
V3 

n = ~2 



VAF , 



(17) 



independent of A, and 



v±_ 



3A 



VAF 



2(A + 1) 

For A >> 1, we note that v± saturates at 

/3 



2^ F - 



(18) 



(19) 



That wj_ saturates with large A is in contrast to the 
behavior of site-centered cases and can lead to rather 
isotropic spin wave cones for the bond-centered case^^ 
despite local microscopic anisotropy. As discussed in the 
next section, for bond-centered stripes with any spacing 
p, vh is independent of A and v± saturates with large A. 
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FIG. 4: Spin wave spectra and intensities for vertical, site- 
centered stripes. All spectra are reported at k y — ir as a 
function of the transverse momentum k x . The frequency uj is 
in units of J a S. Apparent crossings only occur at A = 1 and 
A = 2.5. 



C. Numerical Results 

For most values of the stripe spacing p, the spin wave 
matrices are sufficiently large that one must use numeri- 
cal diagonalizations to obtain the dispersion relations of 
the various modes. From the corresponding cigenfunc- 
tions we can then also calculate the spectral intensity 
(proportional to the dynamic structure factor) that these 
magnon states would contribute to the inelastic neutron 
scattering. Figs. Q] and [S] show the calculated disper- 
sion and scattering intensities for site- and bond-centered 
vertical stripes of various spacings. Our results for site- 
centered stripes are consistent with those of Refi^i. For 
the site-centered case, bands never cross for A < 1. At 
the critical couplings A = 1 and A = 2.5, site-centered 
bands appear to cross. Away from these couplings, verti- 
cal site-centered bands generally repel rather than cross. 
For A = 1, the dispersion is very similar to that of a pure 
antiferromagnet, albeit with different magnetic recipro- 
cal lattice vectors. For any coupling A, as p — > oo, the 
result for a pure 2D antiferromagnet is recovered. For 
p increasing but finite, the number of bands as well as 
the number of reciprocal lattice vectors increases. How- 
ever as p — > oo, all spectral weight is transferred to the 
response of a pure antiferromagnet. 

Fig. [3] shows representative results for vertical bond- 
centered antiphase domain walls with spacings p = 2,3, 
and 4. In this case, the critical point where bands appear 
to touch is at A c ~ 0.56 and is at most only weakly depen- 
dent on p. Away from the critical coupling, bands never 
appear to cross, but rather level repulsion is observed. 
There are other notable differences between the site- and 
bond-centered cases. For one thing, for the same spac- 
ing p, bond-centered configurations yield one more band: 
site-centered configurations have p — 1 bands, whereas 
there are p bands for bond-centered configurations. 

A qualitative difference between the two cases is the 
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(a) Acoustic spin wave velocities for VS3 



FIG. 5: Spin wave spectra and intensities for vertical, bond- 
centered stripes. All spectra are reported at k y = ir as a 
function of the transverse momentum k x . The frequency 10 is 
in units of J a S. 



scaling of the band energies with coupling A. For site- 
centered configurations, all bands increase their energy 
monotonically with the coupling ratio A. This is in con- 
trast with the bond-centered case, where for large 
only the top band is affected by the ferromagnetic cou- 
pling (that is, it increases linearly with A), but all other 
bands saturate as A is increased. The behavior of the top 
band can be understood by considering the spins that are 
ferromagnetically coupled across the domain wall. In the 
top band, these spins precess ir out of phase with each 
other, and the dispersion is dominated by the behavior of 
the effective ferromagnetic dimers, yielding uj — > 2\Jb\S/h 
as \Jb\ — > oo, as shown in Appendix FBI 

An important consequence of the saturation of the 
lower bands as A gets large in the bond-centered cases 
is that the low-energy spin wave velocities alone, v± and 
Uii , cannot readily be used to extract the relation between 
the bare exchange couplings J a and J^. We explore this 
point in more detail in the next section. 

In Fig. H3 we present the spin wave velocities perpen- 
dicular (v±) and parallel (um) to the stripe orientation for 
the acoustic (lowest) bands as functions of the coupling 
constant ratio A. These are compared to the reference 
velocity, vaf of the pure antiferromagnet, which is inde- 
pendent of the coupling A and equivalent to p — > oo. 

While in both the site- and bond-centered cases 
(FigS-UJa) and (b), respectively) the perpendicular veloc- 
ity depends on the coupling ratio, in the bond-centered 
case v± rapidly saturates to a value close to v\i for large 
A. As a consequence, the value of the coupling ratio 
A = J?,/ J a cannot be determined solely by the ratio of 
the acoustic velocities but requires additional informa- 
tion, such as vaf- 

The curves of v± and Vn cross at A = f for the bond- 
centered case, apparently independent of p for the widths 




(b) Acoustic spin wave velocities for VB3 



FIG. 6: Spin wave velocities for (a) VS3 and (b) VB3 as a 
function of the coupling ratio A. The solid curves in panel (a) 
are analytic results for VS3 calculated in Section ff V Bl Sym- 
bols in both (a) and (b) are numerical results. The velocities 
parallel and perpendicular to the stripe direction are equal 
to each other for X — 2/7 and A = 1 in the site- and bond- 
centered cases respectively. Qualitatively similar behavior is 
found for other stripe spacings. 



we have studied. The crossing is at most weakly depen- 
dent on p in the site-centered case, occurring at A = 2/7 
for DS3, and at A = 0.3 for DS4. For all spacings studied, 
we find that in the bond-centered case, v\\ is independent 
of the coupling A and that v±_ rapidly saturates with large 
A. As p gets larger, both of these velocities approach 
vaf- For the VB3 configuration, uii = 0.9v af, indepen- 
dent of A. Notice that the independence of Uii upon A 
and the rapid saturation of v± as A becomes larger than 
1 means that bond-centered configurations can produce 
rather isotropic spin wave conesA* 
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FIG. 7: (Color online) Diagonal site- and bond-centered con- 
figurations, showing even and odd spacing. Dotted lines de- 
note domain walls. Solid parallelograms denote unit cells. 



V. RESULTS FOR DIAGONAL STRIPES 

Fig. depicts representative diagonal configurations, 
for site- and bond-centered domain walls and with even 
and odd spacing. As mentioned in Sec III Bl for a given 
microscopic model, diagonal stripes are more strongly 
coupled across the domain wall than vertical stripes. In 
addition, there are more parameters to consider for site- 
centered diagonal stripes: we must include J c as well as 
Jb (see Fig.0), since both couplings appear to the same 
order if derived from, e.g., a Hubbard- like model. 



ther at (0, 0, w) if the diagonal in-plane stripes lie directly 
on top of each other from plane to plane (meaning there 
is also no net magnetization on a domain wall), or at 
(0,0,0) if the stripes are diagonal within a plane and in 
their correlation from plane to plane. 



B. Analytic results for small p 



As for the case of vertical stripes discussed in Sec. lIVBl 
for small p, it is possible to obtain analytic forms for the 
acoustic dispersion relations for diagonal stripes in both 
the site- and bond-centered cases. 

For the case DS3, the analytic dispersion is 

(^f) 2 /2 = f{k x — k y ) + A 2 f(2 (k x — ky j) 
+ 2 AA C f(k x — k y ) 

+ {X + X c )[f(2k x + k y ) + f{k x + 2k y )] 
+ X[f(3k x )+f(3k y )] , (20) 

(21) 

where A c = | ^f- | , and where the function / is defined as 
f(x) = 1 - cos(x) , (22) 



as in Sec. lTVRl 

The dispersion perpendicular to the stripes, along the 
k = [k x , k x ) direction, is then 



u)(k x , k x 



/ 3k x • 



J S = V8(2A + A C ) |sin(^)|, 
which yields for the velocity in that direction 



3V2A + A c 



(23) 



(24) 



A. Elastic Peak at (0, 0) 

Like their vertical counterparts, bond-centered diago- 
nal stripes can produce new peaks in the elastic response. 
With diagonal stripes the new weight is physically trans- 
parent. For all bond-centered domain walls, nearest- 
neighbor spins are ferromagnetically coupled across the 
wall, and in the diagonal case, nearest neighbor pairs 
along a single domain wall all have their moments point- 
ing in the same direction, leading to a domain wall mag- 
netization. As Fig J7J illustrates, for diagonal stripes 
with evenp, adjacent domain walls have alternating signs 
of the magnetization. But diagonal stripes with odd 
spacing have the same magnetization direction on each 
domain wall. This generically leads to net ferromag- 
netism and a peak at (0,0), unless parameters are fine- 
tuned. In a three dimensional antifcrromagnet, (as may 
happen^ e.g., with weakly coupled planes) domain walls 
are two dimensional (planar), and this peak appears ei- 



which approaches v± — > 2 A vaf for large A, and v± — > 
%VX~ C vaf for large A c . 

In the parallel direction (k x ,—k x ), the dispersion be- 
comes 



1 ( Lo (k x ,-k x ) 

JaS 



= (1 + A + A c + cos k x + X cos (2k x )) 

k 

x (l + 2A + 2Acosfc K )sin 2 y , (25) 

(26) 



which gives 



v/(l + 4A)(2 + 2A + A c ) 



VAF- 



" 2^2 

This approaches Dm — * A vaf for large A, and v\\ 
\/Ac/8 vaf for large A c . 
For the case DB2, the analytic dispersion is 



(27) 



8 




FIG. 8: Dispersion and intensities for DS3 and DS4 along FIG ' 9: Dispersion and intensities for DS3 along (k x ,k x ) 

/, ; \ j- ■,, t i rpi r ••■if (k x ,—2k x ) directions, comparing the effects of Jb and J c . 

(kx,k x ) direction with Jb only. 1 he frequency u is m units of v ' ' ' r b 

J a S. For all plots, J c = . frequency w is in units of J a S. 



(f^) 2 / 2 = 4A(1 + A) + A - AV16(1 + A) 2 +B (28) 
where 

A = (1 — X 2 )f(k x - k y ) (29) 
B = -8(1 + Xff{k x -k v ) 

+ (2-f(k X -ky))[2f(k X -ky) 

- 2f(2k X - 2ky) - f(3k X + ky) ~ f{k X + 'iky)} . 

Perpendicular to the stripes, along k = (k x ,k x ), the ve- 
locity is 



v± = \IJTi VAF ' (30) 

saturating to v± — > Vaf as A >> 1. In the direction 
k = (k x , —k x ), parallel to the stripes, the velocity is 



V]l = ^V^+T v AF . (31) 

C. Numerical Results 

In Fig-El we plot the dispersion and intensities for DS3 
and DS4 along (k x ,k x ) for various values of the coupling 
ratio A = setting J c — 0. (See Fig. ^ for the defini- 
tions of Jb and J c .) Similar results using a J c only model 
(i.e. with Jb = 0) are reported in Refi2^. Our results 
show similar band structures but with critical coupling 
A = 1, which is only half of the J c only model. 

However, the effects of Jb and J c depend upon the 
direction in k space. In Fig. [5] we show the effects of 
varying the couplings Jb and J c for two cuts in momen- 
tum space for DS3. For a cut perpendicular to the stripe 



direction, Jb and J c have more or less the same effect, al- 
though since Jb couples more spins than J c , it has a more 
dramatic effect. Increasing either coupling broadens the 
bandwidth in a roughly linear manner with negligible ef- 
fect on the shape. However, for the cut (k x ,—2k x ), we see 
that the presence of Jb produces inflection points when 
J c = 0, and can produce flat-topped dispersions if J c is 
included as well. 

We show in Fig. ^3 the calculated dispersion relations 
and intensities for the bond-centered diagonal case, for 
spacings p = 2,3, and 4. As in the vertical case, the 
number of bands is equal to p. A striking difference in 
the spectra of odd spacings is seen, as the net ferromag- 
netism in the system changes the low-energy character of 
the spin waves from a linear (antiferromagnetic-like) to 
a quadratic (ferromagnetic-like) dispersion. Rather than 
the band repulsion observed in the vertical case (except 
at finely tuned values of the coupling) , crossing of optical 
bands is generic in the bond-centered diagonal case. Note 
the ability of optical bands to cross, indicating a differ- 
ence in symmetry for the crossing bands. Also evident 
in the dispersion of DBA is the downturn of the acoustic 
band at 2p magnetic reciprocal lattice vectors, twice as 
many as in the odd case. (See Sec. Ill Bl .) This is ex- 
pected because of the doubling of the unit cell necessary 
to accommodate even spacing. Note, however, that spec- 
tral weight is forbidden at these extra reciprocal lattice 
vectors, including the (ir, n) point. 

In Fig. ^2 we plot the spin wave velocities for DS3. 
When Jb and J c are both finite, there is a wide range of 
couplings A for which the spin wave velocities parallel and 
perpendicular to the stripes are nearly equal, while this 
approximate isotropy is confined to a narrow range of A 
if either Jb or J c is zero. Fig. 1121 which presents v± and 
v\\ for the case DB4, shows the characteristic saturation 
of v± with large A for bond-centered stripes. 
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FIG. 11: Velocities parallel (v«) and perpendendicular (v±) 
to the stripe direction, as compared to vaf, for DS3. In 
the first panel, J c = 0, and the velocities are plotted as a 
function of At = .lb I J a. In the second panel, Jb = 0, and the 
velocities are plotted as a function of A c = J c j J a . In the third 
panel, Jb = Jc, and the velocities are plotted as a function of 
A;, = A c . 
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FIG. 12: Velocities parallel (vu) and perpendendicular (v±) 
to the stripe direction, as compared to vaf, for DB4. 



VI. EXPERIMENTAL IMPLICATIONS 

We have shown that for a certain class of nontrivial 
spin orderings on a lattice, the spin wave response is sen- 
sitive to the microscopic placement of the antiphase do- 
main walls. Furthermore, even elastic neutron scattering 
can in principle distinguish site- from bond-centered for 
odd stripe spacings, whether vertical or diagonal. 

While both site- and bond-centered odd width verti- 
cal stripe configurations will produce elastic weight at 
(7r± ^, 7r), only configurations that are phase-shifted from 
the site-centered configuration (e.g. a bond-centered con- 
figuration) are capable of producing weight at (0, n), and 
the observation of this peak along with peaks at (ir±j, it) 
would rule out a site-centered vertical configuration. A 
similar ferromagnetic peak-ie., at (0,0)-would rule out 
site-centered diagonal stripesiii 

Figs. and 1121 illustrate another important implica- 
tion for experiments: The transverse spin wave velocity 
fx in the acoustic band saturates for large A in the bond- 
centered case for both vertical and diagonal stripes. In 
fact, all but the top band in the bond-centered case satu- 
rate and become independent of J& for large Jb- As noted 
above, this unfortunately means that an estimate of 
A = J b / J a cannot necessarily be discerned directly from 
the ratio i?xA>|| but requires either independent knowl- 
edge of whether the stripes are site- or bond-centered, or 
an appropriate estimate of the bare coupling (from, e.g. 
vaf). 

A prominent piece of phenomenology in the cuprates 
is the the "resonance peak" observed in neutron 
scattering which is the presence of extra scatter- 
ing weight appearing at (n, tt) at finite frequency, typ- 
ically of order AOmeV. One proposal is that this may 
be due to spin waves crossingi^Si We note that for ver- 
tical stripes, spin waves generically repel and appear to 
cross only at finely tuned values of the coupling. For 
site-centered configurations, this corresponds to A = 1 
and A « 2.5, while for bond-centered configurations, 
the critical coupling is near A ~ 0.56 . However, a fi- 
nite energy resolution measurement would not be able to 
distinguish actual crossings from near crossings. In the 
bond-centered case with large A, the first optical mode 
has more weight at (71, 7r) than the acoustic band, which 
would tend to leave the weight near the elastic incommen- 
surate peaks (7r ± jj , 7r) rather disconnected from what 
might be called a "resonance peak" in this configuration. 
We also note that our calculations show that band cross- 
ings are more generic in the presence of diagonal stripes 
than vertical stripes. 

The nickelate compound Lai.69Sro.3iNi04 shows evi- 
dence from neutron scattering of diagonal stripes with 
spacing p — 3££ As Sr is substituted for La, holes are 
doped into the MO2 planes. Neutron scattering has 
been used to map out the acoustic spin wave dispersion 
for this material. The data reveal rather isotropic spin 
wave cones, i.e. that v±_ and v\\ are rather similar, with 
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v ± w (1.03±0.06)uaf andwn w (0.86 ±0.06)t)AFr where 
w J 4i7' is the acoustic spin wave velocity of the undoped an- 
tiferromagnet. For the DS3 state, if we include only Jb 
or only J c , we find no coupling strength A for which these 
two relations can be simultaneously satisifed. The pres- 
ence of the two couplings together, as shown in Fig. 1771 
with Jb = J c , can account for the proper relationship 
among the velocities, but only for a small range of rather 
small coupling ratio. As a general trend, we find ap- 
proximate isotropy of the spin wave cones to be more 
robust for bond-centered stripes (in both vertical and di- 
agonal cases), and so one might suspect bond-centered 
stripes could be responsible for the near isotropy of the 
spin waves in this material. However, as we have shown, 
the DB3 configuration yields a ferromagnetic spin wave 
dispersion, which is certainly not supported by the data. 

In the related compound La2Ni04.i33i£, signatures of 
spin stripes have been detected in neutron scattering. 
"Incommensurate peaks" are observed to persist up to 
a temperature T m , above which magnetic peaks indica- 
tive of stripe structure can be regained by application 
of a 6T magnetic field. The field-induced stripe spacing 
(both above and slightly below T m ) is smaller than the 
zero-field stripe spacing observed below T m . As noted 
by the authors^ the ferrimagnetic response is naturally 
explained by bond-centered stripes. In the high tem- 
perature field-induced stripe phase, the diagonal stripes 
have spacing p = 3. Our results in Fig. 1101 suggest that 
this field-induced transition should be accompanied by a 
dramatic change in the low-energy spin wave dispersion, 
from linear to quadratic. 

We have also shown that (as in the site-centered case24) 
the number of bands in a bond-centered configuration is 
set by the number of spins in the unit cell, rather than 
by the spacing p. Generally, for both vertical and diag- 
onal stripes, site-centered stripes have {p — l) spin wave 
bands, and bond-centered stripes have p bands. The ex- 
ception is the case of diagonal site-centered stripes with 
odd spacing p, which has \{p — 1) bands. An experi- 
mental consequence of this is that for a given value of p, 
bond-centered stripes have p spin wave bands, whereas 
site-centered stripes have at most p — 1 bands. Although 
not yet observed experimentally, this means that the 
upper bands can also be used to distinguish site- from 
bond-centered stripes. Finding p bands along with in- 
commensurate peaks indicative of spacing p would rule 
out site-centered stripes. For diagonal odd width stripes, 
the threshold is even lower. For, e.g., DS3, only one 
spin wave band is expected, whereas for DB3, we expect 
to find three bands. The observance of a second band 
(or eqivalently a spin wave crossing) for diagonal p = 3 
sripes would rule out a site-centered configuration. Of 
course, negative evidence is dicier, and the observance of 
the smaller number of bands cannot distinguish the two, 
as it cannot rule out the possibility that the top band is 
too faint to be observed. 



VII. CONCLUSIONS 

In conclusion, we have studied regular arrays of an- 
tiphase domain walls in two-dimensional Heisenberg an- 
tiferromagnets and find that their location relative to 
the lattice-i.e., whether they are site-centered or bond- 
centered-produces distinct effects which may be mea- 
surable in a diffraction probe such as neutron scatter- 
ing. In particular, arrays of odd-width, bond-centered 
antiphase domain walls generically produce more elas- 
tic peaks than site-centered stripes. In addition, bond- 
centered stripes generically produce more bands than 
site-centered stripes. We further find that low-energy 
spin wave velocities are not always directly related to the 
exchange couplings in the model, and in particular for 
bond-centered configurations, rather isotropic spin wave 
cones are predicted for a wide range of parameters. 
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APPENDIX A: SPIN- WAVE METHODS 

We rewrite the Hamiltonian Eqn. using the ladder 
operators: 

H = \ E JrAS z r S z r , + \(StS- + S-S+)]. (Al) 

<r,r'> 

We now replace the spin operators by Holstein-Primakoff 
(HP) bosons^ 

S+ = V2Sai 
Sr = V2Sa+ 

Sf= S- aj a, (A2) 
for odd sites i occupied by a spin up, and 
S+= V2Sa+ 

S? = -S + a+a, (A3) 

for even sites i occupied by a spin down. Here, i labels 
each spin within a unit cell, i.e. i — 1, 2, • • • , N, where 
N is the number of spins in the unit cell. We use odd i 
to represent S z =T spins and even i for S z =[ spins. We 
Fourier transform the bosonic operators via 

a t (k) = — V a r e ikr , 

r£odd i 

V " r- 

r^even % 
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Finally, we get the Hamiltonian in momentum space 

H = ^A lja +(k) aj (k) 

ij 

+ 5E[^at(k)a+(k)+S*.a i (k)a i (k)], (A5) 

ij 

where A = A + and B T = B. 

The quadratic Hamiltonian l|A5|l can be diagonized via 
a canonical symplectic transformation 33 T, b = Ta, using 
the bosonic metric 

where / is the N x N indentity matrix. This leads to 
H (k) - ]T K (kK (k)6 Q (k) + \u> a (k)] . (A7) 

a 

We now consider the structure factor. Only S x and S a 
contribute to the inelastic part of the structure factor. In 
terms of HP bosons, 



(A8) 



sjl( ^K+(-k) + ai (k)] 

i£odd 

+ J2 [a+(k)+a 4 (-k)] 



We then substitute Eqn. I|A8(I into the structure factor 
and keep only the creation operators {bf (k), b£ (k), • • •}, 
which connect the vacuum to singly excited states. This 
gives 

S m (k,u; a ) = 2j2\<f\SZ\0>\ 2 S("-"f) 
f 

= s\<i\(£ <*ib+)\o>\ 3 

i 

= S\J2"i\ 2 > (A9) 



where ot% is the ith component of the (orthonormalized) 
eigenvector \a > of the Hamiltonian using the bosonic 
metric, corresponding to eigenvalue u> a . 



APPENDIX B: DIMERIZED SPIN MODEL 



We consider an isolated system of two spins with ferro- 
magnetic coupling J;,. In the ground state, the two spins 
are aligned. When the spins tilt a bit, each produces an 
effective field acting on the other. Using the classical spin 
wave method, we have 



*- 4* *ft <■»> 

Ignoring the change in S z , the x, y components of the 
two spins satisfy 



dSf 
~dt~ 
dS\ 
dt 



dS£ 

dt 
dS£ 

dt 



JbS 



(Sf 



sf) 



(B2) 



Integrating yields Sf' v = —S%' v + c, where c is a constant 
of integration. Since we allow only S z to have a constant 
component, c = 0. Taking the second derivative of Sf, 
we find 



d 2 Sf 



J b S dSf dSj x 



dt 



dt 



_JbSy j Q(QX Q x \ 

Ti 



.J 2 h s 2 



[Sf - Sf ] 

QX 



(B3) 



which is a harmonic oscillator equation. If we set 
Sf (x, t) = u(x)e luJt , we see that the oscillation frequency 
is 



ft 



which recovers the large J& limit of Eqn. El 



(B4) 
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